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Abstract 

With the discovery of evidence for neutrino mass, a vivid gamma ray sky 
at muIti-TeV energies, and cosmic ray particles with unexpectedly high energies, 
astroparticle physics currently runs through an era of rapid progress and moving 
frontiers. The non-vanishing neutrino mass establishes one smooth component 
of dark matter which does not, however, supply a critical mass to the Universe. 
Other dark matter particles are likely to be very massive and should produce high- 



> 

QQ , energy gamma rays, neutrinos, and protons in annihilations or decays. The search 

for exotic relics with new gamma ray telescopes, extensive air shower arrays, and 
underwater/-ice neutrino telescopes is a fascinating challenge, but requires to un- 
^\ , derstand the astrophysical background radiations at high energies. Among the 

0^ ' high-energy sources in the Universe, radio-loud active galactic nuclei seem to be 

,S^ • the most powerful accounting for at least a sizable fraction of the extragalactic 

f^, gamma ray flux. They could also supply the bulk of the observed cosmic rays at 

O ' ultrahigh energies and produce interesting event rates in neutrino telescopes aim- 

■4-J ' ing at the kubic kilometer scale such as AMANDA and ANTARES. It is proposed 

C^ . that the extragalactic neutrino beam can be used to search for tan lepton appear- 

ance thus allowing for a proof of the neutrino oscillation hypothesis. Furthermore, 
a new method for probing the era of star formation at high redshifts using gamma 
rays is presented which requires new-generation gamma ray telescopes operating 



5^ , in the 10-100 GeV regime such as MAGIC and GLAST. 



1 Introduction and practical definition of high-energy 
astroparticle physics 

Central to modern astronomy is the dark matter problem and it is commonly 
believed that its solution will trigger major advances in particle physics and cos- 
mology 1^]. So far dark matter is only known through its gravitational effects, 
but the understanding of the nature and origin of dark matter requires to obtain 
more direct information about its mass and interactions. Cosmology and particle 
physics qualify weakly interacting massive particles (WIMPs) with masses between 
100 GeV and a few TeV as likely candidates. The WIMPs violently annihilate with 
their anti-particles in rare collisions or they could be unstable. These modes lead 
to secondary gamma rays and neutrinos which can be detected on Earth H. The 



suspected large WIMP mass then corresponds to gamma ray and neutrino energies 
in excess of 100 GeV. In addition to the dark matter particles, there may be other 
relics from the early Universe, such as quintessence, vacuum energy, or topolog- 
ical defects. Quintessence is a slowly rolling scalar field connected with massive 
bosons, and this could also lead high-energy phenomena, although no worked-out 
model exists to my knowledge. The recent discovery of accelerating expansion using 
SNIa as a tracer of cosmic geometry seems to make a strong case for quintessence 
[P, 0]. Similarly, if the scalar field has settled to some false (meta-stable) vacuum, 
the energy density of this vacuum could drive the deceleration parameter away 
from n/2. Topological defects preserve false vacuum in pointlike (monopoles), 
one-dimensional (strings), two-dimensional (domain walls), or higher dimensional 
space-time structures. They are topologically stable, but have a variety of ways to 
communicate to our world in addition to their gravitation. E.g., they can dissipate 
into GUT-scale bosons (~ 10^^ GeV) which are unstable themselves and fragment 
into jets consisting of gamma rays, neutrinos, and protons at ultra-high energies 
p|. After their propagation through intergalactic space, electromagnetic cascading 
and secondary particle production shift most energy injected by these exotic pro- 
cesses to much lower energies where the energy release competes with that due to 
ordinary astrophysical sources. 

In order to identify new physics phenomena, it is therefore of crucial impor- 
tance to obtain a com,plete inventory of the astrophysical high-energy sources which 
act as a background for these searches. This defines the high-energy astroparticle 
physics program from a practical point of view and follows the logic inherent to the 
general astronomical exploration of the sky to cover the entire range of wavelengths 
with comparable sensitivity. 

Among the non-thermal sources in the Universe, radio-loud active galactic nu- 
clei (AGN) seem to be the most important energetically. There are other interest- 
ing sources, such as gamma ray bursts (GRB) and clusters of galaxies, but their 
non-thermal energy release does not come close that of AGN. There are some in- 
triguing complications arising through calorimetric effects, since the intracluster 
medium in clusters of galaxies surrounding AGN confines escaping relativistic par- 
ticles for some time and thereby gives rise to a secondary luminosity tied to the 
energy release of the AGN and the cooling time scale of the intracluster medium. 
The radio-loud AGN come in various disguises, depending on the orientation of 
their radio jet axes and the properties of the circum-nuclear matter in their host 
galaxies. The most extreme versions are radio galaxies with the radio jet axes 
almost in the plane of the sky and the blazars with the radio jets pointing close to 
the line of sight to the observer which leads to a dramatic flux increase owing to 
special relativistic effects (the so-called Doppler boosting). In Sect. 2 it is argued 
that radio-loud AGN can be expected to produce the entire extragalactic gamma 
ray background from an energetical point of view. Owing to beaming statistics, 
the number of unresolved sources responsible for most of this background should 
not be too large. Indeed the flux from the ~ 50 resolved sources in the fiux-limited 
EGRET sample already equals a sizable fraction (of the order of 15%) of the extra- 
galactic background fiux. With the next-generation gamma ray telescopes MAGIC 



1^ (giant 17ni air-Cerenkov telescope) and GLAST |^ (space-borne silicon strip 
detector) it will be possible to probe deeper into the astrophysical source popu- 
lation producing the extragalactic background below 100 GeV thereby narrowing 
the range of opportunity for particle physics models of exotic processes producing 
gamma rays. In fact, the gamma ray background below 100 GeV provides a good 
measure of the entire electromagnetic energy release during the history of the Uni- 
verse, since gamma rays from remote sources cascade down to the energy range 
below 100 GeV which is shown in Sect. 3. It is pointed out in Sect. 4 that a MAGIC 
observing campaign for high-redshift gamma ray sources can be used to probe the 
era of star formation and the evolution of the optical-ultraviolet metagalactic ra- 
diation field back to redshifts of ~ 5. 

The recent discovery of multi-TeV emission from Mrk 421 and Mrk 501 H, the 
measurement of their spectra using the HEGRA air-Cerenkov imaging telescopes 
P], and the improved measurements of the extragalactic infrared background |10|, 



make a strong point in favor of accelerated protons in extragalactic radio sources 
which is shown in Sect. 5. The following Sect. 6 discusses the immediate implications 
that the radio-loud AGN could well produce the observed cosmic rays at highest 
energies and high-energy muon neutrinos. Finally, in Sect. 7 it is pointed out that 
the extragalactic muon neutrino beam is likely mixed with tau neutrinos fO] which 
leads to very interesting experimental signatures, such as the disappearance of the 
Earth shadowing effect at ultra-high energies and the appearance of tau leptons in 
underwater /-ice detectors. 

2 Origins of extragalactic background radiation 

Inspection of Fig. 1 shows an interesting pattern in the present-day energy density 
of the diffuse isotropic background radiation consisting of a sequence of bumps 
each with a strength that is decreasing with photon energy. The microwave bump 
is recognized as the signature of the big bang at the time of decoupling with its 
energy density given by the Stefan-Boltzmann law usk = aT'^. The bump in the 
far- infrared is due to star formation in early galaxies, since part of the stellar 
light, which is visible as the bump at visible wavelengths, is reprocessed by dust 
obscuring the star-forming regions. The energy density of the two bumps can be 
inferred from the present-day heavy element abundances. Heavy elements have a 
mass fraction Z = 0.03 of the total mass density p.^, and were produced in early 
bursts of star formation at redshift Zf by nucleosynthesis with radiative efficiency 
e = 0.007 yielding 



1 + -Zf 

Inserting plausible parameter values one obtains 



(1) 




nn. ~ 6 X 10-3 tlll^ i— ^ eV cm-3 (2) 



for the sum of the far-infrared and optical bumps. Probably all galaxies (except 
dwarfs) contain supermassive black holes in their centers which are actively accret- 
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Fig.l: Sketch of the present-day energy density of the extragalactic 
radiation background from radio waves to gamma rays. The dashed 
hne shows the expected AGN contribution to the low-energy diffuse 
background from the average quasar spectral energy distribution. 

ing over a fraction of tagn/i* ~ 10^^ of their lifetime implying that the electromag- 
netic radiation released by the accreting black holes amounts to 



"agn 



n, 



s -- 1.4 X lO"*^ eV cm" 



(3) 



adopting the accretion efficiency eaccr = 0.1 and the black hole mass fraction 
Mbh/-M^* = 0.005 [|l^]. Most of the accretion power emerges in the ultraviolet where 
the diffuse background is unobservable owing to photoelectric absorption by the 
neutral component of the interstellar medium. However, a fraction of Ux/ubh ~ 20% 
taken from the average quasar spectral energy distribution |jl^ shows up in hard X- 
rays due to coronal emission from the accretion disk to produce the diffuse isotropic 
X-ray background bump with 

2.8 X 10"^ eV cm"^ 



u^ 



(4) 







Jets with non-thermal 7-ray emission show up only in the radio-loud fraction 
^ri ~ 20% of all AGN and their kinetic power roughly equals the accretion power 
|15|. Hence one obtains for the background energy density due to extragalactic jets 



0.2 



— ] 2.8 X 10"^ eV cm~3 
0.2/ 



(5) 



If unresolved extragalactic jets are responsible for the diffuse gamma ray back- 
ground, this requires a particle acceleration efficiency given by 



u 



U\ 



?rad'Wj 



(6) 



where Uj denotes the total (kinetic + magnetic + randomized relativistic particle) 
energy density in extragalactic jets. Inserting the energy density of the observed 
extragalactic gamma ray background]^ one obtains a limit for the acceleration effi- 
ciency 

Cacc > o.ise-;^ (ll) ' (7) 

which is of the same order of magnitude as the 13% efficiency required for super- 
nova remnants to produce the Galactic cosmic rays. Accelerated protons achieve 
this high radiative efficiency, if they reach energies of up to 10* TeV. In the next 
section it is shown that protons at such high energies cannot go unnoticed, they 
produce interesting gamma ray spectra owing to the photo-production of secon- 
daries. Some of the protons turn into neutrons due to vr"*" production and can leave 
the jets without adiabatic losses. These particles would have just the right flux 
to produce the observed extragalactic cosmic rays dominating the local spectrum 
above lO^*'^ eV, as will be shown in Sect. 6. 

3 Cascading and gamma ray calorimetry 

Gamma rays of energy E can interact with low-energy photons of energy e from the 
diffuse isotropic background over cosmological distance scales I producing electron- 
positron pairs 7 -|- 7 — s- e"*" -|- e~ , if their energy exceeds the threshold energy 

2(mec^)^ ^fl + zy^f E Y' 



where // denotes the cosine of the scattering angle |lq|. The 7-ray attenuation e~'^ 
due to pair production becomes important if the mean free path A becomes smaller 
than I, i.e. if the optical depth across the line of sight through a sizable fraction of 
the Hubble radius obeys r = l/X > 1. For the computation of r one first needs to 
know the pair production cross section 



^(1-.^) 



2/3(/32-2) + (3-/?^)ln(i±|) 



(9) 



where /3 = ^1 — 1/7^ with 7^ = e/eth, and where ctt denotes the Thomson cross 
section |jl^. Then one needs the geodesic radial displacement function dl/dz = 
-jj-[(l -|- z)E{zy\^'^ to compute the line integral from z = to some z = Zo- For a 
cosmological model with Q. = 1 and A = the function E{z) simplifies to (H-z)^'^. 
Hence one obtains the optical depth 

<'^° , dl /•+! , 1 



T^^{E,Zo)= / dz—- / djj. — -— / den^{e){l + z) ajy{E,e,^,z) 

Jo dz J-l I Jeth 

= — - / dz^i + zf'"^ dx- denb{e)a.,j{E,e,x-l,z) (10) 

Ho Jo Jo 2 Je . 



^ Note that the flux in the gamma ray background observed by CGRO is close to the bolometric 
gamma ray flux of the Universe, since pair attenuation and cascading must lead to a steepening 
of the background spectrum above 20 — 40 GeV jl6|, [l7[ . 
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Fig. 2: The diffuse isotropic microwave-to- ultraviolet background. 
Solid curve: 10th order polynomial interpolation of observational 



data (|20|, g2|, |21|1, and references in |16|). 



adopting a non-evolving present-day background density nb, i.e. n'^{z, e')de' = (1 + 
z)^n\y{e)de where the dash indicates comoving-frame quantities. The simplifying 
assumption that the photon density transforms geometrically corresponds to the 
situation in which an initial short burst of star formation at Zf > Zq produced most 
of the diffuse infrared-to-ultraviolet background radiation. This simple assumption 
is replaced by a more realistic one in Sect. 4. Fig. 2 shows the spectrum of the 
low-energy diffuse background used to solve Eq.(lO) numerically. 

Figure 3 shows the resulting t{E, z) = 1 (omitting the subscript hereafter) curve 
for the microwave-to- ultraviolet diffuse background spectrum shown in Fig. 2. It is 
obvious that 7-rays above ~ 10 — 50 GeV cannot reach us from beyond redshifts 
of z = Zf = 2 — 4. Higher energy 7-rays can reach us only from sources at lower 
redshifts (e.g. 7-rays with energies up to 10 TeV have been observed from Mrk 501 
at z = 0.033 in accord with Fig.3 ^). 

Corollary I: // the extragalactic gamma ray background originates from unre- 
solved sources distributed in redshift similar to galaxies, its spectrum must steepen 
above ~ 30 GeV due to 'y-ray pair attenuation. 



Here is has been tacitly assumed that the 7-rays which have turned into electron- 
positron pairs do not show up again. This is, in fact, not quite true, since the pairs 
are subject to inverse-Compton scattering off the microwave background thereby 
replenishing 7-rays. The 2.7 K background is more important as a target than the 
shorter wavelength background, since there is no threshold condition for Thomson 
scattering contrary to pair production and since 2.7 K photons greatly outnumber 
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Fig. 3: The 7-ray horizon t(E, z) = 1 for the low-energy background 
spectrum shown in Fig. 2. Cosmological parameters are h = 0.6, 
= 1, and Qa = 0. For a general discussion of pair attenuation, see 



reference [23|. 



the latter. The inverse-Compton scattered microwave photons turn into 7-rays of 
energy 

'1 + z' 



Eu 



10 



SOGeV; 



MeV 



(11) 



conserving the energy of the absorbed 7-ray which corresponds to a constant 
E'^dN/dE, i.e. the expected slope of the differential spectrum is about -2 (-2.1 
observed). A small amount of energy is lost to lower frequency synchrotron emis- 
sion, if magnetic fields are present in the interagalactic medium. 

Corollary II: Energy conservation in the reprocessing of 'j-rays from higher 
to lower energies by pair production and subsequent inverse-Compton scattering 
produces an approximate dN/dE oc E^^ power law extragalactic gamma ray back- 
ground between ^ 10 MeV and ~ 30 GeV. 



Evolution of the met agalactia optical- ultraviolet ra- 
diation field 



4.1 A simple model based on the observed "effective" star forma- 
tion rate 

Consider an effective^ cosmic star formation history p*(2;) denoting the production 
rate per unit volume of mass which has formed to stars at a redshift of z. Such star 
formation histories have been inferred from galaxy counts in the Hubble Deep Field 
[p4| . Since the present-day infrared background is strong enough to aborb gamma 
rays in the TeV range in the local Universe far from the peak of the star formation 
history, its evolution in the past is rather irrelevant in this context. However, the 
present-day optical radiation background scaled back to the peak of the star forma- 
tion history at a redshift of Zb ~ 1.5 implies gamma ray attenuation in the 20 GeV 
regime from sources at this redshift. Over this distance scale the evolution of the 
background becomes important, since it is gradually produced by the forming stars. 

For the gamma ray attenuation only the evolution of the number density n(e, z)d€ 
of the background photons is important which relates to the photon production rate 
h oc p/e in the following way: 

n(6, z)de = £' dz'Me', z')^ (^) ' de' (12) 

The evolving background must be normalized to yield the observed present-day 
radiation background 

f^i fif Q 

n{e,0)de= / dz''ii{e',z')^ {I + z'y^ de' (13) 

Jo dz 

As a simple example consider a burst of star formation at a high redshift h oc 
6{z — Z{). Inserting this in Eq.(12) and combining with Eq.(13) we obtain 

n{€,0)de = {l + z)-^n{e',z)de' (z < z^) (14) 

which represents a constant co-moving density background density where the 
(1 -|- z)^ term reflects the geometric scaling of the cosmic volume. To obtain a 



realistic parametrization of n{z,e')de' we approximate the Madau curve |24] as a 
broken power law 

h{z, e')(l + z')-^ oc ("Mf) ") oc (1 + z)°-i (15) 

with a = om = 3.8 for < z < 1.5 = Zb and a = /3m = —4.0 for Zb = 1-5 < 
z < W = Zf. We also investigate a star formation rate which exhibits a plateau 
beyond Zh- 

Equation (1) enters the formula for the gamma ray optical depth: 

rz di r+l roo 

t{Eo,z)= / dz—- / dfi^l — fi) / den{e,z)a{E,e, ^) (16) 

Jo dz J_l Je^i^ 



^The term "effective" means that only the star formation rate inferred from photons which have 
made it through possible obscuring dust clouds are of relevance for the build-up of a metagalactic 
radiation field. 



Note that cosmology enters through dl/dz (which depends on 0,17^, and Ho), 
not through n(e, z) for a given parametrization oi p^,{z). However, the parametriza- 
tions p^:{z) must also satisfy observational constraints such as number counts and 
the present-day diffuse background which themselves depend on cosmology. Turn- 
ing this around it means that one must find the cosmology parameters for which a 
measured gamma ray horizon (i.e., the curve t{Eo,z) = 1) and the star formation 
history data come into mutual consistency. 

The gamma ray horizon from Eq.(16) is shown in Fig. 4 using the low-energy 
background spectrum template shown in Fig. 2. The template was used to normal- 
ize the evolving background such that is identical to the template at z = and 
scales to higher redshifts acording to Eq.(12). The fact that there is a redshift 
with a maximum star formation is very important. If power law evolution of the 
background emissivity were to continue all the way in the past, one could easily 
infer power law solutions for the scaling of n(e, z) which are more shallow than 
(1 + z)"^ as in ref. |lf:], but such solutions become unrealistic beyond zy,. 



4.2 Extragalactic gamma ray background 

The origin of the observed diffuse isotropic gamma ray background is unknown. 
The spectral shape and flux density suggest that unresolved faint radio loud AGN 
are responsible for this background, similar to the situation in the X-ray band 
where deep observations have revealed that faint AGN are responsible for more 
than 90% of the background emission. The EGRET-type radio-loud AGN seem 



contribute not more than ~ 25% j25| to the extragalactic gamma ray background. 
The uncertainties about the faint end of the gamma ray luminosity function in the 
EGRET band allow for a larger contribution from the general class of radio-loud 
AGN. According to beaming statistics, the flux-limited EGRET sample of AGN 
is dominated by highly beamed sources with a rather fiat luminosity function. A 
much fainter, less beamed population with a steeper luminosity function is likely 
to fill in the remaining 75%, at least this seems very plausible considering the 
energetics of radio jets as was shown in Sect. 2. I strongly expect that the flat- 
spectrum/steep-spectrum classes are mirrored in different populations of gamma 
ray sources. The nearest steep-spectrum radio sources would have been detected 
by EGRET even if they were faint (with their gamma ray luminosity roughly a 
factor of ~ 50 larger than the 5 GHz luminosity). However, with a lower compact- 
ness for instrinsic gamma ray absorption oc L/R the steep-spectrum radio sources 
could emit most of their gamma ray power above the EGRET range. It is up to 
new air-Cerenkov telescopes with threshold energies above 10 GeV and GLAST to 
probe this proposal. 

If the extragalactic jets are indeed responsible for most of the gamma ray back- 
ground, it is straightforward to investigate the effect of pair attenuation on the 
spectrum of the background. The precise shape of the spectra of the individual 
sources at high redshifts is rather unimportant owing to the effects of cascading 
discussed in Sect. 3. Adopting a power law gamma ray spectrum per source with 



The cosmological gamma ray horizon 

f2=1 ; f2^=0; h=0.5; z^=^ .5; z,=1 0, a=3.8; P=-4.0 (CDM) 
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Fig. 4: Gamma ray horizon due to interactions with an evolving 
metagalactic radiation field as computed from the effective star for- 
mation rate. Sources below the horizon curve suffer no significant 
pair- attenuation along the line of sight. The grey band indicates the 
uncertainty of the horizon as estimated from the range of interpola- 
tions allowed between observational upper and lower bounds of the 
flux of the present-day optical-UV diffuse background. Note that 
the metagalactic radiation field before the maximum of the cosmic 
star formation rate (indicated by the dashed line) is too weak to 
significantly attenuate gamma rays below 20-40 GeV resulting in a 
near-constant optical depth. The light solid line shows the effect of 
a star formation rate with an extended plateau which causes the op- 
tical depth to continue to grow with redshift beyond z = 1.5. The 
horizontal lines indicate the effective threshold energies for various 
air-Cerenkov telescopes. It is emphasized that triggering below 20- 
40 GeV, which can be achieved by the MAGIC telescope, is crucial 
for probing the star formation era. Such an investigation is comple- 
mentary to studies of galaxies at high redshifts, since it is addition- 
ally sensitive to diffuse sources of optical-UV photons such as would 
be arising from exotic particle decays (one possible scenario for the 
reionization epoch). 
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the average slope of the resolved EGRET sources 

extending from Ei = 10 MeV to £^2 = 1 TeV and taking into account the luminosity 
density evolution ^(z) of AGN, we obtain a good approximation of the present-day 
background energy density 

u{E) = —EIe (18) 

c 

from the equation 

u{E) oc dz^^{z){l + z)-*B ^ J ' e ^^C[E,z] (19) 

Jo dz \ El J 

where dt/dz is given by 

'^^''^'^ " (1 + z)./ri(T+ zf + {l-VL- VLk){1 + zY +11 ^^°^ 

and the function C\E^ z] for the effect of pair attenuation can be approximated as 

C[E,z]=e~^ (21) 

with Et{z) denoting the solution of the equation t{E,z) = 1 (the gamma ray 
horizon). The result is shown in Fig. 5. It will be possible with GLAST to find 
whether the diffuse gamma ray background indeed turns over in this shallow fashion 
or continues as a power law into the 100 GeV domain. In the latter case, the gamma 
ray background would have to be due to some local source population [26|. 



5 Comparison of proton blazar predictions with ob- 
served multi-TeV spectra 

In the previous sections arguments based on energetics have been used to favor 
extragalactic jets as the sources of the gamma ray background. This requires that 
the jets radiate a sizable fraction of their kinetic energy in gamma rays when inte- 
grated over their lifetimes. Since the cooling of relativistic particles increases with 
their energy, a high radiative efficiency is equivalent with high energies. For elec- 
trons, Lorentz factors required for a high radiative efficiency are at least 7e ~ 10^, 
and for protons 7p ~ 10^/(1 + 500u^/uB)^e ~ 10^°. The Lorentz factor for pro- 
tons may seem outrageously high, but in a statistical acceleration process such as 
Fermi acceleration with the balance between energy gains and losses determining 
the maximum energy, such high energies are an inevitable consequence of the accel- 
eration theory. Moreover, particles with energies in excess of 10^^ eV are observed 
in the local spectrum of cosmic rays. Their energy is too high for the gyrating 
particles to be isotropized in the Galactic disk, so that an extragalactic origin is 
very likely. The energy requirements can be converted to an energy supply rate for 
extragalactic sources, and this requires sources as strong as radio galaxies. Thus, 
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Fig. 5: The gamma ray horizon enforces a shallow turnover beyond ~ 
30 GeV in the spectrum of the extragalactic gamma ray background 
which is only weakly dependent on cosmology (results from further 
investigations by T. Kneiske will soon be reported elsewhere.). 

if it is not the radio sources themselves, some other source must be able to pro- 
duce a relativistic proton distribution with an enormous energy flux reaching these 
ultrahigh-energies. It has been suggested that GRBs may do that, but there are 
strong arguments against that proposal presented in Sect. 6. Here we concentrate 
on the consequences of accelerating baryons in radio jets to these extremely high 
energies where radiative losses become important. A few years ago, a model coined 
the proton blazar model has been presented in which accelerated protons have been 



assumed in addition to the accelerated electrons |27]. This model made definite 



predictions about the multi-TeV spectra of nearby blazars [28| which can now be 
compared with observations of Mrk 501 and Mrk 421 obtained with the HEGRA 
air-Cerenkov telescopes. In fact, this was the only model prior to the observations 
which made any quantitative prediction of multi-TeV emission from these sources. 

An important, although not necessary, assumption of the original model is that 
the magnetic field pressure energy density is in equipartition with that in relativistic 
particles implying that synchrotron cooling dominates over Compton cooling (since 
the photon energy density remains below that in particles). This affects accelerated 
electrons as well as secondary electrons at ultrahigh-energies. Evaluating a simple 
conical jet geometry shows that typical blazars are optically thin to gamma rays 
up to the TeV range. Unsaturated synchrotron cascades initiated by accelerated 
protons interacting with the synchrotron photons from the accelerated electrons 
are computed as the stationary solution of a coupled set of kinetic equations which 
is then Doppler boosted to an appropriate observer's frame. A series solution is 
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found employing Banach's fixed point theorem which can be physically interpreted 
as a series of superimposed cascade generations. The cascade generation of gamma 
rays emerging in the TeV range on the optically thin side has a spectral index 
s ~ 1.7 (differential spectrum Iq) steepening by q = 0.5 — 0.7 above TeV. The 
index a is the energy index of the optical synchrotron photons which act as a 
target for both gamma rays and protons. The reason for the break is the onset of 
intrinsic pair attenuation characterized by the escape probability I^ = PqscIo for a 
homogeneously mixed absorber and emitter 

The shape of the multi-TeV spectrum is therefore not sensitive to changes in the 
maximum energy and can remain constant under large-amplitude changes of the 
flux associated with changes in the maximum energy. 

The observed multi-TeV spectrum is modified by the quasi-exponential pair 
attenuation due to collisions of the gamma rays with photons from the infrared 
background. A recent evaluation of this background based on direct measurements 
obtained from COBE data in the far-infrared, and inferred as a lower limit from 
number counts based on ISO observations of the HDF shows that this effect com- 
pensates the shallow downward curvature discovered by the HEGRA collaboration 
in the spectrum of Mrk 501 Q. If the gamma rays were due to inverse-Compton 
scattering, the shallow curvature is difficult to understand for a number of rea- 



sons given in ref. [22|. The most important of them is that one must expect the 
accelerated electrons not to be able to reach energies much higher than 10 TeV. 
An inverse-Compton spectrum produced by these electrons would therefore have 
to show significant curvature approaching this maximum energy adding to the in- 
evitable curvature due to gamma ray interactions with the infrared background 
photons. There are some rumours that quantum gravity effects could possibly sup- 
press pair production over intergalactic distances, but that remains highly specula- 
tive. I consider the agreement between the proton blazar prediction and observation 
very promising for the model, albeit minor discrepancies must be expected, since 
the proton blazar model is highly simplified in order to avoid too many free pa- 
rameters and to be predictive. Therefore I take the freedom to speculate about 
the emissions associated with the gamma rays, viz. cosmic rays and high-energy 
neutrinos in the following sections. 

6 Neutrino and cosmic ray predictions 

The photo-production of pions leads to the emission of neutrons and neutrinos. 
The neutrons decay to protons, and such extragalactic cosmic rays suffer energy 



losses traversing the microwave background [29|. At an observed energy of 10^^ eV, 
the energy-loss distance is Ap '^ 1 Gpc owing to pair production. This distance 
corresponds to a redshift Zp determined by Ap = (c/Ho) Jq'^ dz/[{l + z)E{z)] where 

E{z) = [0(1 + zf + nR{l + zf + n^l ^ with [7 -F [7/j -h J^A = 1- Almost indepen- 
dent on cosmology, the resulting value for Zp is given by Zp = /i5o/(6 — /150) — O.2/150 
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Fig. 6: Comparison of predicted and observed flux density spectrum 
in the multi-TeV range for Mrk 501. The thin solid line shows the 



spectrum published in |28]. The upper thick solid lines show this 
spectrum scaled to the 300 GeV flux levels during the two obser- 
vation epochs where the air-Cerenkov data indicated by the solid 
and open symbols were obtained with the HEGRA telescopes. The 
dashed line shows the spectrum without the effect of the assumed 
marginal interagalactic gamma ray attenuation due to interactions 
of the gamma rays with metagalactic infrared radiation. More re- 
cent HEGRA observations with higher statistical significance show 
some downward curvature in the 10 TeV range which may be at- 
tributed to a stronger attenuation which is in line with new analyses 



of COBE data and ISO galaxy counts in the Hubble Deep Field [10|. 
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where /150 = Ho/50 km s^^ Mpc~^. Therefore, when computing the contribution 
of extragalactic sources to the observed cosmic ray flux above 10^^ eV, only sources 
with z < Zp must be considered. Assuming further that extragalactic sources of 
cosmic rays and neutrinos are homogeneously distributed with a monochromatic 
luminosity density ^(z) oc (1 + z)^'^'' where A; ~ 3 for AGN |3^, their contribution 
to the energy density of a present-day diffuse isotropic background is given by 

^ ^ io cdz Ho Jo {l + z)^E{z) ^ ' 

where 2^ = 2 denotes the redshift of maximum luminosity density. The factor 
(1 + z)^'^ accounts for the expansion of space and the redshift of energy. For 
a simple analytical estimate of the effect of energy losses on the proton energy 
density at 10^^ eV, we collect only protons from sources out to the horizon redshift 
Zp ~ 0.2 for 10^^ eV protons, whereas neutrinos are collected from sources out 
to the redshift of their maximum luminosity density Zm. This yields the energy 
density ratio for neutrinos at an observed energy of ~ 5 10^^ eV and protons at 
10^^ eV 

^.(0) ^ a^-{l + zf-VE{z)dz _ 

np(0) £^{l + zf-VE{z)dz ^ ^ ^ 

using ^ ~ 0.3 from decay and interaction kinematics, and considering an open 
Universe with E{z) = (1 + 2;) and a closed one with E{z) = (1 + z)'^'^. 

Fig. 6 shows exact energy-dependent results for Vt = 1 from a full Monte-Carlo 
simulation employing the matrix doubling method of Protheroe &: Johnson |^] and 
using the model A neutrino spectrum from the original work [33|. The associated 
gamma ray flux corresponds to the observed background flux above 100 MeVJj. The 
neutrino flux is consistent with the bound given in ref. [32], although it is possible 
to have extragalactic neutrino sources of higher neutrinos fluxes without violating 
the observed cosmic ray data as a bound [^4|. Note that there are a few cosmic 
ray events at energies above 10^'^ eV which are difficult to reconcile with an origin 
in extragalactic radio sources, since the radio galaxies are typically at such large 
distances that pion production quenches their spectrum above 10^^'^ eV. However, 
cosmic ray particles from the few closest radio galaxies deflected by magnetic fields 
could possibly explain these events. If not, they might originate from the decay of 
still higher energy particles, such as the gauge bosons produced at cosmic strings 
[^ indicating new physics. 

7 Neutrino oscillations and event rates 

The neutrino flux shown in Fig. 7 corresponds to a very low muon event rate even 
in a km^ detector which is of the order of 1 event per year and per steradian. 
This event rate could be increased if there is additional neutrino production due 
to pp-interactions of escaping nucleons diffusing through the host galaxies which is 

^ A recent paper by Waxman and Bahcall 132] refers to the neutrino flux from model B in the 
original work which was given only to demonstrate that hadronic jets cannot produce a diffuse 
gamma ray background with an MeV bump (as measured by Apollo and which is now known to 
be absent from a COMPTEL analysis) without over-producing cosmic rays at highest energies. 
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Fig. 7: Comparison of proton (solid line) and neutrino fluxes (dotted 
lines, from top to bottom v 



fij i^fij 



and fe) from the proton blazar 
model (Monte-Carlo computations and figure by R.J. Protheroe). 
The open symbols represent the observed cosmic ray flux. 



difficult to predict due to their unknown magnetic fields and turbulence level. The 
reason for the low rate is that the neutrino spectrum is extremely hard, a differential 
proton spectrum of index Sp = 2 photo-producing pions in a synchrotron photon 
target also with differential index Ssyn = 2 yields a differential neutrino spectrum 
of index 3,^ = 1 {dN/dE ex E^^) up to some very high energy. The spectrum may 
be more shallow, if the target photons have a spectral index of 0.7 as suggested for 
Mrk 501, thereby increasing the number of lower energy neutrinos while keeping 
the bolometric flux the same. Because of the long lever arm from 10^'^ TeV to 
1 TeV, a factor of ~ 100 increase in the event rate would result from this effect. 

At this point the discovery of neutrino mass announced by the Super-Kamiokande 
collaboration |0] comes in changing the situation in a major way. A deficit of at- 
mospheric muon neutrinos was observed with Super-Kamiokande at large zenith 
angles with the most likely explanation being a full-amplitude oscillation of muon 
flavor eigenstates to tauon flavor eigenstates across the Earth at GeV energies. 
While this would make long-baseline experiments searching for the appearance of 
the tauon in a muon neutrino beam with laboratory beams extremely difficult (if 
not impossible), it qualifies the expected extragalactic sources of muon neutrinos 
as an ideal neutrino beam. The energies are high enough to produce tauons on the 
mass-shell and the distance large enough to obtain full mixing. Since tauons de- 
cay before interacting in the Earth and since the Earth is opaque to tau neutrinos 
above ~ 100 TeV, a fully mixed extragalactic muon neutrino beam must initiate 
tauon cascades in the Earth shifting the tauon neutrino flux down to energies of 



100 TeV 1 35] and obliterating the Earth-shadowing effect [36| that makes the 
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niuon solid ange very narrow at high energies. 

The neutrino oscihations would have another important consequence. Current 
data suggested a maximal mixing between muon and tauon flavors and a mass 
difference given by 

Am% = ml - m\ = 5 x lO'^ eV^ (25) 

The mass difference between electron and muon flavored neutrinos inferred from 
the solar neutrino deficit is orders of magnitude less, and this implies that the 
neutrino masses could be highly degenerate if they are at the eV level. Using the 
limit rriy^ < 5 eV, the maximum allowed combined neutrino mass would be 

m^ ^ 3m^^ < 15 eV. (26) 

Inserting this into the Cowsik-McClelland bound one obtains the maximum con- 
tribution to the (hot) dark matter of the Universe 

O, < i^/j-2 ~ 0.4 f-^V' . (27) 

91.5 V0.65y ^ ^ 

Due to free streaming of the neutrinos at the time of recombination, density fluc- 
tuations of this dark matter component would be wiped out on scales less than 

A^ = 30(V0.65)"^ Mpc (28) 

and would therefore have no consequence for the dark matter inferred from studies 
of galaxy halos or galaxy clusters which yield O = 0.3ib0.1 [^]. Structure formation 
simulations exclude il^ to be larger than 0.15 |^^ which is in agreement with the 
above upper limit and both evidences together rule out neutrinos as the dark matter 
which could supply a critical mass to the Universe. 

8 Discussion and summary 

The paper highlights in a personally biased way on a few developments in high- 
energy astroparticle physics, rather than to give a review of all the activities in the 
field which encompass a much wider scope from the origin of cosmic rays to quan- 
tum gravity and which involve a truly impressive number of experimental efforts. 
From my point of view, the identification of one component of dark matter using 
the discovery of the atmospheric neutrino anomaly with Super-Kamiokande repre- 
sents a major achievement, and I have highlighted its consequence that neutrinos 
cannot close the Universe. Furthermore, this discovery is not sufficient to prove 
the neutrino oscillation hypothesis conclusively. An experiment is needed which 
shows the appearance of the tau lepton in a beam of muon neutrinos, and I have 
shown that neutrinos due to proton acceleration in extragalactic sources would 
be ideally suited as a beam for a tau-appearance experiment in one of the major 
neutrino telescopes which are under construction. From an astrophysical point of 
view, the oscillations are also important for another reason, since they remove the 
Earth shadowing effect which suppresses the response of a neutrino telescope to 
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extraterrestrial sources of very high-energy neutrinos. This makes the discovery of 
neutrinos from radio-loud AGN much more probable in a few years. 

Other kinds of dark matter such as right-handed neutrinos, supersymmetric par- 
ticles, magnetic monopoles, or other topological defects and their associated gauge 
bosons, are likely to be very massive. It is a formidable task to find other than grav- 
itational evidence for this dark matter due to the electromagnetic (photons), weak 
(neutrinos), and strong (protons) couplings of its decay and annihilation products. 
Inevitably this task requires to understand all the astrophysical background radi- 
ations at high energies. The astrophysics governing this energy domain itself has 
proven to be a fascinating realm. A good example is the discovery of multi-TeV 
emission from nearby blazars using the HEGRA imaging air-Cerenkov telescopes. 
The observational findings frustrate the worldwide elite in theoretical astrophysics 
and seem to provide key insights into the extraordinary physics driven by weakly 
accreting black holes. 

If the recent determinations of the infrared background from ISO galaxy counts 
(lower limits) and COBE/DIRBE and FIRAS are correct, then the steepening 
in the multi-TeV spectrum of Mrk 501 observed with the HEGRA air-Cerenkov 
telescopes is due to the predicted gamma ray attenuation in collisions with these 
infrared background photons. Predictions of the attenuation process for sources at 
higher redshifts depend strongly on the evolution of the metagalactic optical-UV 
radiation field. It will therefore soon be possible to independently probe the star 
formation history using 10-100 GeV gamma rays from extragalactic high-redshift 
sources when lower threshold gamma ray telescopes such as MAGIC and GLAST 
are available. The method is sensitive to truly diffuse photons between the galaxies 
and, when compared with calculations based on the observed star formation rate 
in early galaxies, allows to test the hypothesis whether the background photons 
originate from the stars or other sources (possibly connected with the reionization 
epoch). 

The multi-TeV spectra from nearby blazars prediced on the basis of the proton 
blazar model are in accord with the observations if the effect of pair attenuation due 
to the extragalactic infrared background is taken into account. This is surprinsing, 
since the model gives stationary spectra, whereas the observed flux is highly vari- 
able. Nevertheless, the spectra seem to remain rather constant in the multi-TeV 
domain. There are several properties which are not explained in the framework 
of a stationary quasi-homogeneous model by contraction, such as the apparently 
different variability pattern of Mrk 421 in the GeV and TeV ranges. This does not 
argue against the proton acceleration hypothesis, since one could construct inhomo- 
geneous, non-stationary models to explain this phenomenology. The same is true 
for the problem of short-term variability, which may require to describe the pas- 
sage of shocks traveling through narrowly-spaced inhomogeneities. It is amazing, 
however, that a general feature expected from models based on electron acceler- 
ation is certainly not seen in the data of Mrk 501, viz. the change of the upper 
cutoff energy with varying flux. This can only mean that the cutoff is indeed due 
to intergalactic attenuation and that the electron maximum energy is much higher 
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than 20 TeV which is difficult to understand in the presence of synchrotron and 
inverse-Compton losses. Another peculiar finding is that the multi-TeV emission 
in Mrk 501 is accompanied by 100 keV synchrotron photons in some epoches, but 
not in all of them. Moreover, in Mrk 421 100 keV synchrotron flares have not been 
observed in spite of flaring TeV emission. So the story is not simple and likely to 
continue, brainstorming is always desired when confronted with a new phenomenon. 

An interesting corollary from hadronic models of extragalactic gamma ray 
sources is that they would also emit cosmic rays and neutrinos at about equal 
luminosities (within a factor of a few) . If the cosmic ray flux emitted by hadronic 
accelerators equals that of the observed cosmic rays above 10^^ eV, the associated 
gamma ray power from these sources is enough to produce the observed extra- 
galactic gamma ray background above about 100 MeV. The gamma ray power is 
larger than that in cosmic rays, since the cosmic rays lose energy traversing the 
low-energy background radiation fields and most sources have high redshifts. It 



has been proposed in ref. [39| that GRBs are responsible for the highest energy 
cosmic rays, but more recent GRB observations indicate that most of them have 
high redshifts which is expected if they trace star formation. The strong evolution 
of the GRB luminosity density would then rule out GRBs as possible sources of 
the highest energy cosmic rays, since their cumulative gamma ray flux is far be- 
low the putatively extragalactic gamma ray flux. Although the muon event rate 
in neutrino telescopes from hadronic extragalactic gamma ray sources supplying 
the highest energy cosmic rays is low, neutrino oscillations lead to tau cascades 
canceling the Earth shadowing effect thereby increasing the detection probability. 
Tau lepton appearance in the neutrino telescopes would constitute the first direct 
measurement of the tauon (which so far has only been inferred from momentum 
conservation) and would prove the neutrino oscillation hypothesis. 
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